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Prom random polynomials to symplectic geometry 

Steve Zelditch 



Abstract. We review some recent results on random polynomials and their generalizations in complex 
and symplectic geometry. The main theme is the universality of statistics of zeros and critical points of 
(generalized) polynomials of degree N on length scales of order -7= (complex case), rcsp. ^ (real case). 



1. Introduction 

This is a short survey of some results of P. Blehcr, J. Neuheisel, B. Shiffman and the author on random 
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Motivation 



to study random polynomials and their generalizations in geometry comes from several sources: 

• Classical Analysis: Value distribution theory of polynomials and analytic functions is a classical 
topic. Computable examples may exhibit non-generic patterns of zeros (or other values) and one 
would like to understand the typical distribution. One forms ensembles of analytic functions by 
defining the coefficients to be independent random variables with a gi ven d i strib ution. One can then 
study expected behaviour, almost sure behaviour and so on (see e.g. [jKacj, |lo|, |(|). 

• PDE: Spherical harmonics of degree N are examples of eigenfunctions of the Laplacian on a compact 
Riemannian manifold. They are restrictions to the sphere S m of homogeneous harmonic polynomials 
of degree N on R m+1 . One would like to know about nodal lines, critical points, sup norms (etc.) of 
general Laplace eigenfunctions. Studying features of random spherical harmonics gives insight into 
the 'typical' properties of eigenfunctions and avoids pathologies such as occur in U [JN). Analo gucs 
of spherical harmonics of degree N on general compact Riemannian manifolds are linear combinations 
of eigenfunctions Aip\ = A 2 tp\ with A 6 [JVlog N, (N + l) log(iV + l)]. References on random spherical 
harmonics include |Be| , [n| , [v|, ^el ); for random combinations of eigenfunctions, sec [Ze2|. 

• Quantum Chaos: Eigenfunctions of quantum chaotic Hamiltonians are well modeled by random 
polynomials in regard to distribution of zeros or critical points, to sizes (sup- norms or L p -norms), 
to quantum expectation values, and in other respects (see, among others, |ABST| , |BBL| , [HKZ|. 
NoVo, Zel, Ze2|). This is analogous to the similarity between eigenvalues of random matrices and 
eigenvalues of quantum chaotic systems. 

• Algebraic Geometry: Holomorphic sections s £ H°(M, L N ) of the Nth power of an ample line bundle 
L — > M over a Kahler manifold (M, oj) are quite analogous to homogeneous polynomials of degree N, 
and coincide with such polynomials when M = CP m , L — 0(1) ( CP™ = complex projective m-space, 
0(1) is the hyperplane bundle (cf. JGH| )). The simultaneous zero set Z S1 ,... iSk of k holomorphic 
sections defines a codimension k algeraic submanifold of M; one would like to know the 'almost sure' 
properties of such a submanifold. 

• Symplectic Geometry: Almost- holomorphic sections s £ Hj(M,L N ) of ample line bundles L — » M 
over almost-complex symplectic submanifolds (M, J, u>) in the sense of [Don, BoGu] are very similar 
to holomorphic sections in the complex case. Under (difficult) transversality conditions, they have 
applications in symplectic geometry analogous to those in algebraic geometry. 
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In this article, we will refer to the setting of holomorphic (or almost-holomorphic) sections as the 'complex 
case', and the setting of eigenfunctions of Laplacians as the 'real case'. The complex wave functions live 
on phase space while the real eigenfunctions live on configuration space. The main theme of our work has 
been the universality of statistics of zeros and critical points of polynomials of large degree N on small 
length scales (of order in the complex case, resp. in the real case). We have only considered compact 
manifolds, and only Gaussian or spherical measures on their spaces of polynomials. Two universality classes 
have emerged: (i) the 'Heisenberg class' in the complex case (connected to the Heisenberg group), and (ii) 
the 'Euclidean class' in real case (connected to the Euclidean motion group). It should be mentioned that 
the 'real case' has many other meanings in the literature on random polynomials (e.g. real polynomials and 
their real zeros). 



2. Mathematical Tools 



We give a quick summary of some basic tools and methods that are used below. We will concentrate on 
general ideas and refer to [fehZe , 5hZe2 , BSZ3 for detailed expositions. 



2.1. Vector spaces of large dimension. In defining our ensembles of polynomials, we will be dealing 
with a sequence (Hn, (, )n) of Hubert spaces of increasing dimension rfjy = dimHN, where N is the 'degree of 
the polynomial', a large integral (semiclassical) parameter. The dimension is given by a (Hubert) polynomial 
of the form ~ a^N" 1 in the complex case (with to = dimcM) and by a function of polynomial growth 
dN ~ aoN rn ^ 1 in the real case (with to = dirngM). 

In our applications, the spaces H.n will be one of the following. 

2.1.1. Spherical harmonics and real eigenfunctions. We denote by A the standard Laplacian on S m 
and by TtN(S m ) the space of spherical harmonics of degree N on S m . They are the eigenfunctions of A of 

/ N+m-1 

eigenvalue \ m ,N = N(N+m—l) and form a real vector space of dimension d m .N = 2 n+"i-i 

y m — 1 

More generally, we may consider the Laplacian A of any compact Ricmannian manifold (M, g) . In place of 
spherical harmonics of degree N, we partition the spectrum of a/A into intervals [iVlogTV, (N + l) log(iV+ 1)] 



(the reason for the longer length is given in [Ze2|), and let Hn denote the span of the eigenfunctions with 
eigenvalues in the iVth interval. Linear combinations of such eigenfunctions are of course not eigenfunctions of 
A, but they behave like polynomials of degree N and their random linear combinations give a replacement for 
random spherical harmonics. Their use for modelling quantum ergodic and quantum mixing eigenfunctions 



is discussed [Ze2, HKZ 



2.1.2. Holomorphic sections of positive line bundles. For any Kahler manifold (M, u>) of complex dimen- 
sion m there exists a holomorphic hermitian line bundle (L, h) — > (M, oS) whose Ricci curvature Ric(h) = lu. 
L is called positive since it posseses a metric of positive curvature, i.e. oj(X, JY) defines a Riemannian 
metric. We denote by L N the iVth power of L and by H°(M,L N ) the space of holomorphic sections. Its 
dimension — dimH°(M, L N ) is given by the Hubert polynomial dpj — Cl N m + • • • for sufficiently 
large N, where • • • represent the lower order terms. We equip M with the volume form dV u = ^-r and 



H°(A1, L N ) with the inner product ||s|| 2 = J M h(s(z), s(z))dV. For background we refer to [GH 
In the simplest case of Riemann surfaces, examples include: 

• M = CP 1 ,^ = 0(1), L 2 = TCP 1 ,/! = hps,0J — lufs (Fubini study hermitian metric and curvature 
(l,l)-form). That (CP 1 , ujfs) is positively curved in the usual Riemannian sense is equivalent to posi- 
tivity of TCP 1 . 7?°(CP 1 , 0(2N)) may be interpreted as the space of holomorphic vector fields of type 
(-§^) N ■ More simply put, sections are homogeneous holomorphic polynomials s = X^=o c k^o z i °f 
degree N in two complex variables. Such polynomials are known as the S'C/(2)-ensemble. 

• M = H 2 /r,L = T*M,hps — hyperbolic metric. Hyperbolic surfaces are negatively curved in the 
Riemannian sense, so their tangent bundles are negatively curved, and their co-tangent bundles are 
positively curved. Holomorphic sections are holomorphic differentials of type (dz) N . 

• M = C/Z 2 ,wo = dz A dz. The complex torus is flat in the Riemannian sense, so neither its tangent 
nor cotangent bundles are positively curved. The 'quantizing line bundle' with curvature loq is rather 
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the bundle O whose sections are the classical theta-functions. The sections of Q N are known as 

theta-functions of level N. 
2.1.3. Almost holomorphic sections. Symplectic almost-complex manifolds (M,J,u>) possess a similar 
but analytically more complicated geometric quantization as spaces Hj(M,L N ) of 'almost- holomorphic' 
sections. They are defined by a D complex over the S^-bundle X due to Boutet de Monvel -Guillemin 
[BoGu]. What we need to know about these spaces is that their orthogonal projectors II have the same 
scaling asymptotics as in the complex case, if one works in suitable (Heisenberg) coordinates | ShZe2| | . 

2.2. Gaussian measures and spherical measures. We will restrict attention to two related ensem- 
bles: 

• Gaussian ensembles: We fix an orthonormal basis {fj} of TIn and write functions as orthonormal 
sums / = YljZi c jfj- We then define the (complex) Gaussian measure by 

7 = e" ll/l|2 X>/, i.e. 7 = e~ |c|S dc. 

More generally we fix a symmetric matrix A on C rfjv with positive (semi-)definite imaginary part and 
define: 

e -(A- 1 c, c ) 

(1) 7A= (2^detA dc ' 

Gaussian ensembles come in both real and complex flavors. In the real case, the exponents acquire 
factors of 1/2 and the denominator acquires a square root. 

• Spherical ensembles: We denote by STLn = {/ S TLn ■ ||/|| = !}■ We then equip STLm with the 
uniform (Haar) probability measure v^q . 

We will denote the expected value of a random variable X with respect to the Gaussian ensemble by 
E 7 (X) (resp. E„(X) for the spherical ensemble). 

These two ensembles are equivalent in the sense that the two large dimension limits give equivalent 
results when scaled properly. A more precise formulation goes as follows: 

Let Tn : R dN — ► R k , N — 1,2,..., be a sequence of linear maps, where oIm — > oo. Suppose that 
j^T N T N -> A. Then T N *v dN -> j A . 

2.3. Sequences of random polynomials. We are often interested in sequences of polynomials {sjv} 
chosen independently and at random from TCn from cither a Gaussian or spherical ensemble. We thercform 
form the product probability space (HocMoo), defined by 

Hoc = Hi x H2 x • • • x Hn x • • • , jUoo = x^ =1 ^i N , (li n = 7at or u N ). 

When we say that a sequence of polynomials {siv} of increasing degerees does something almost surely, we 
mean that the set of such sequences has measure one in this product ensemble. 

2.4. Szego kernels. Our results depend on the fact that certain statistical properties of polynomials 
can be expressed in terms of the reproducing kernels 11^ (x, y) (orthogonal projections) of the Hilbert spaces 
TLn- They are known as Szego kernels, and are essentially the same as the 'coherent states' of the physics 
literature. The local structure of the Szego kernel is given by the following scaling asymptotics: 

Theorem 2.1. As N — > oo, we have: 



BSZ2, ShZe2] Complex case: 



rWzn + — — z n + — — ) - J-e^-^e^-sCN'+H 2 )/! + _L Pl ( uv . Zn ) + ...\ 



The leading order term is the Szego kernel for the reduced Heisenberg group, whence the name 'Heisen- 
berg class. ' 
Real case: 

n-^u n ( Xq + ^,.x + ~ r(^)(J^)^ J^(\u -«!) + ■.., 

where J v is the Bessel function of the first kind of order v (whence the name 'Euclidean class'.) 
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The proof of the scaling asymptotics in the complex holomorphic case [BSZ2] is based on the Boutet- 
de-Monvel- Sjostrand parametrix for the Szego kernel, which is valid for positive line bundles. A similar 
paramctrix was constructed in the symplectic almost-complex case [3hZe2], and the scaling asymptotics 
were derived from it. In the real case of S m , the scaling asymptotics are closely related to the 'Mehler- Heine 
formula'. The terms 'Heisenberg class' and 'Euclidean class' suggest infinite dimensional Gaussian ensembles 
related to representations of the Heisenberg and Euclidean motion groups. 



3. Distribution of zeros and critical points 

We now state some results on the zeros and critical points of random generalized polynomials. Let 
( Sl ,...,s k ) € H°(M, L N ) k or let (s u ...,s k ) £ H k N in the real case. When k = 1 we omit the subscript. We 
will use the following notation. 

• Z s — {x : s(x) = 0} denotes the zero set of s, and \Z S \ denotes the m — fc-submanifold Riemannian 
volume density induced by u> in the complex case or by g in the real case. We further denote by 
\\Z S \\ the mass of \Z S \ and define the probability measures Z s = rrg4r ■ When one takes m sections 
(or functions) in dimension m, then the simultaneous zeros are almost surely a discrete set and the 
measure Z s is the normalized sum of delta-functions at the zeros. 

• C s — {z : Vs(z) = 0} denotes the critical point set of s. In the holomorphic case, V is the holomorphic 
connection compatible with h. We note that C s is almost surely a discrete set. We define |C/| = 

• By the density of zeros at degree TV we mean the coefficient K^ k (z) of the measure K^ k (z)dV = E|Z S |, 
i.e. f M ipE\Z s \ = E J M ip\Z s \ for tp G C[M). Similarly, we denote by K^ nt ' (z) the density (relative 
to the given volume form) of EC S . 

• More generally, we define the pair correlation densities of zeros (resp. critical points) by K^' k (z 1 , z 2 )dV 

E(\Z S \ x \Z S \), resp. K^ k ' N {z 1 ,z 2 )dV = E(\C S \ x |C S |). They are densities of measures on M x M. 
Roughly speaking, the two-point correlation gives the probability density of finding a pair of zeros (or 
critical points) at {z l ,z 2 ). More generally, there are n-point correlation functions, but for simplicity 
we only consider n = 1,2. 

3.1. Statement of results. We have results on several levels: expected values, almost sure behaviour, 
and scaling asymptotics. The following theorems are valid on any complex or almost-complex symplectic 
manifold (M, a>), equipped with a hermitian complex line bundle of curvature ui. 

Theorem 3.1. In the complex case, the density of zeros, resp. critical points, satisfies: 
ShZej] Complex zeros: K* k (z)dV = tu k + 0(±). 



ShZe3] Complex critical points: There exists a universal constant 7 m depending only on the dimen- 
sion such that Kl rlt ' N {z)dV = r y m ^N m + 0(N m ^ 1 ), In particular, the expected number of critical 
points is given by E#C S = -i m Vol L0 {M)N m + 0{N m ~ 1 ). Here, VoL(M) = Cl( ^, )m is the volume of 
(M,lo). The density of critical points is therefore universal. 

Thus, zeros and critical points tend to concentrate in regions of high curvature. Similar results should 
hold in the real case. In the case of S m , such density results are obvious since they must be rotationally 
invariant. 

3.2. Almost sure distribution of zeros. A deeper question on distribution of zeros is whether se- 
quences of individual sections tend to have equidistributed zeros. 

Theorem 3.2. We have: 



ShZefl For a random sequence in the complex codimension k case, {sat}, Z Sn — * u k almost surely. 

] For a random sequence of spherical harmonics on S m ,m > 6, we have Z SN — > dvol almost surely. 
The same is true in (Cesaro) mean for dimensions < 6. 
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3.3. Universality and scaling of correlations. The next level of results concerns the statistics of 
zeros and critical points on the length scale -j= (complex case) or ^ (real case) . For brevity we only describe 
the results in the complex case. Upon magnifying a small ball B_d_ (zq) around an arbitrary point zq of this 

radius, one loses track of the specifics of the geometrical setting and obtains universal limiting correlations. 
More precisely, such a universal limit occurs if one chooses the coordinates properly. 

Theorem 3.3. Scaling limits of correlations of zeros and critical points are universal. That is: 

1 k n ( -> (V 1 7 2 ) 



N 2k ^VVlv'Viv 

in the sense of measures. 

The limit correlation function denoted K^^z 1 , z 2 ) is unique within a universality class for each statistic 
under consideration. In the complex and almost-complex cases, there exists a limit correlation function 
K 2 ^ ros '° (z 1 , z 2 ) of zeros, and another limit correlation function K^^ 8 ' 00 (z 1 , z 2 ) for critical p oints. T here 
are analogous results in the real case. The limit correlations are explicitly computable. In |BSZ1 | and 



elsewhere we give explicit formulae for low values of m and graph the results. The result for zeros in 
dimension one agrees with the formula of Hanny |Han | on CP 1 , as it must since the result is universal. 



We now briefly describe the elements of the proof. The details differ, but the principles are the same, for 
complex and real cases, and for zeros or critical points. Hence we concentrate on zeros in the complex case. 

3.3.1. Step One: Relating correlations and joint probability distributions. Following an idea due originally 
to Kac and Rice in the case of real polynomials of one variable, we express the correlation measures in terms 
of the joint probability distribution of the random variables s(z 1 ), . . . , s(z n ), Vs(z 1 ), . . . , Vs(z"). This JPD 
is defined by 

Df :=D%(x,Z,z)dxdZ = (J x )*u If , M*) = (s{z),Vs(z)) ) 
i.e. it is the push- forward of the Gaussian measure under the linear jet map J(z). 

The desired expression for correlations in terms of the JPD is given by the following generalization of 



the Kac- Rice formula Kac| to the geometric setting of this article: 



Theorem 3.4. jBSZlj , |BSZ3| , $hZe2H In the case of correlations of complex zeros, we have: 



« n 

Kg(z) = J dtD»(o,z,z)i[det(een 



P =i 

Analogous formulae exist for critical points, but involve the jet maps (Vs(z), V 2 s(z)). The real case 
is similar to the complex case, but is somewhat more complicated because TIjv is oscillatory rather than 
exponentially decaying Jn| . 

3.3.2. Step two: scaling asymptotics of the JPD. The JPD is a (generalized) Gaussian measure on the 

complex vector space of 1-jets: = Ja n (z) > where the covariance matrix A N (z) = g J 1S gi ven 

in terms of the Szego kernel and its covariant derivatives, as follows: 

A =W) = ^n w (^,0;^' ) 0), B={Bl, ql ) = £vjn*(z*, 0;**\0), 

C=(C™,) = ±Vl%Il N (zPAz p ',0), P,p' = l,...,n, q,q> = l,...,2m. 

Here, Vj, respectively V 2 , denotes the differential operator on X x X given by applying V 9 to the first, 
respectively second, factor. The link between the JPD and the Szego kernel stems ultimately from the fact 
that n^v is the covariance matrix of 7 on TIn- Using the scaling asymptotics of the Szego kernel we obtain 
(in the complex case): 

Theorem 3.5. ( HShZe2[ , Theorem 5.4) With the same notations and assumptions, we have: 



'(z 1 z ") * s a un ^ versa ^ Gaussian measure, and A (z/yN) — » A a °(z) 
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is given in terms of the Szego kernel for the Heisenberg group: 

, . ml ( A°°{z) B°°{z)' 



0), 



c x {L) m \ C°°(z) 



U^( Z p,0;zp' ,0) for 1 < q < m 

for rn + 1 < q < 2m 

-zP)(z^,-z^))Ilf(z p ,0;z p ',0) for l<q,q'<m 

for max(<7, q') > m + 1 

The scaling limit thus gives the correlations between zeros in the 'Heisenberg ensemble', an infinite 
dimensional Gaussian ensemble. In the real case (i.e. S m ), the limit correlations coincide with those in an 
ensemble related to the Euclidean motion group, and involving the Bessel kernel JnJ. 

3.4. Hole probabilities. Another application of the scaling applications is to 'hole probabilities', i.e. 
the probability that a ball B t (zq) of radius r around a point zq G M is zero-free. The following result 
combines our scaling asymptotics and Sodin's reformulation (and substantial simplification) of Offord's large 
deviations results on hole probabilities for entire analytic functions in the plane |So| . It is based on the 
Poincare-Lelong formula, so at this time of writing it has only been proved for one holomorphic section, and 
no comparable results have been proved for random spherical harmonics. 

Theorem 3.6. Let P N (D;z ) = Prob{s G H°(M,L N ) : Z s n B_^(z ) = 0}. Then there exists positive 

constants C\,C<i such that, for any N, Pjf(D;z ) < C\e~ CriD . 

4. Quantum ergodicity and random waves 

We end with a brief discussion of the connection between random polynomials and quantum chaos. 
The intutitive idea that quantum chaotic eigenfunctions should resemble 'Gaussian random waves' seems to 
have been first suggested by M. V. Berry |b[| . A precise formulation of this random wave model and some 



numerical results are given in |ABST] |HR 



Random spherical harmonics on S m , or random combinations of eigenfunctions of Laplacians on gen- 
eral compact Riemannian manifolds as described above, provide a rather different random wave model for 
quantum chaotic eigenfunctions. To motivate the model, we recall that the diagonal sums of squares 

(3) S p (X) = S2 \(A<pj, ipj) - o~ A \ p , (with cT A = / <J A d[i), 

,,\ • A J S'M 

and their off-diagonal analogues, are used to characterize eigenfunctions as quantum ergodic, quantum mixing 
and so on. Here, A is an observable (zeroth order pseudodifferential operator), a a is its principal symbol. 
For quantizations of ergodic systems, S P (X) — > as A — > oo, and it is of some interest to measure the rate 
and to relate it to the class ical dynamic s. In the random spherical harmonics model one has the following 
rate (for similar results see |Ze2 , ShZe|): 



Theorem 4.1. (|Zel|, Lemma (2.15)) Let {<fNj} be a random orthonormal basis of spherical harmonics 



ofS\ Then: E{S 2 {N)) = M^k?) Js-S* \°T (0 " <^l 2 ^(0] + 0(1/N 2 ). 



The optimist may conjecture a similar rate for quantum chaotic eigenfunctions. In | Zel | it is further 
proved that almost all orthonormal bases of spherical harmonics are quantum ergodic, and this was improved 
to quantum unique ergodicity by VanderKam |jv[| . On the level of quantum mixing, however, random 
spherical harmonics do not provide a good model, but random combinations of Laplace eigefunctions on 



generic Riemannian manifolds do; that was the motivation for studying the model in [Ze2|. Quantum 
mixing systems involve off-diagonal sums like (^) with constant gaps between eigenvalues. Conversely, if the 
eigenfunctions satisfy (||) and the analogous off-diagonal estimates for all gaps, then the system is classically 
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mixing. Hence the random wave model is not a good model for ergodic systems which fail to be mixing (cf. 



[ HKZ| ). 

A further relation between Gaussian random waves and quantum chaotic eigenfunctions does not seem to 
have been explored, even numerically. Let {^a} be eigenfunctions of a quantum chaotic system, and consider 
the local rescaling ifi\(xo + j). The rescaled eigenfunction is an eigenfunction of the rescaled operator, which 
is asymptotically Euclidean. Hence tp\(xo + j) can be expanded asymptotically in terms of plane waves, 
and one might ask how the frequencies are distributed. 
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